Compressibility of the Two-Dimensional infinite-U Hubbard Model 
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We study the interactions between the coherent quasiparticles and the incoherent Mott-Hubbard 
excitations and their effects on the low energy properties in the U = oo Hubbard model. Within the 
framework of a systematic large-N expansion, these effects first occur in the next to leading order in 
1/N. We calculate the scattering phase shift and the free energy, and determine the quasiparticle 
weight Z, mass renormalization, and the compressibility, ft is found that the compressibility is 
strongly renormalized and diverges at a critical doping S c = 0.07 ± 0.01. We discuss the nature of 
this zero-temperature phase transition and its connection to phase separation and superconductivity. 
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In recent years, there has been growing interests in 
the physics of doped Mott-insulators in connection with 
high-T c copper-oxide superconductors. In the absence of 
a natural small parameter, the relevant models of strong 
correlation have been extended and studied under large 
symmetry groups (large N) or large dimensions (large 
d). A generic feature of strong correlation is the coexis- 
tence of coherent quasiparticles 0] and the broad inco- 
herent Mott-Hubbard excitations [ 3 that carry the main 
part of the spectral weight at small doping. It has been 
shown in the t-J model that the systematic large-N ex- 
pansion in the slave boson formalism provides a transpar- 
ent non-perturbative description of both the low-energy 
Fermi-liquid like quasiparticles (^] already present in the 
large-N limit, and the incoherent Mott-Hubbard features 
at next-to-leading-order in 1/N [Q]. 

In this paper, we study corrections to the low energy 
properties due to the effects of the interactions between 
quasiparticles and the incoherent Mott-Hubbard excita- 
tions by a complete calculation of the free energy and the 
single particle Green's function to next-to-leading order 
in 1 /N. This has not been understood properly because 
of the difficulty involved in calculating the corrections 
to the mean-field parameters. For simplicity, we shall 
consider the U — oo Hubbard model with the spin sym- 
metry group generalized from SU(2) to SU(N), although 
the physics discussed here pertains to models that include 
superexchange interactions such as the t-J model. This 
model has been solved for N = oo. The ground state is 
a Fermi liquid at finite hole concentrations and exhibits 
a Brinkman-Rice transition at half-filling ||. We find 
that the interactions represented by the 1 /N fluctuations 
are very strong near half-filling, giving rise to a divergent 
compressibility at a finite critical doping 5 C — 0.07 ±0.01 
below which the Fermi liquid phase becomes unstable. In 
contrast to the Brinkman-Rice transition at half-filling 
in the large-N limit, the quasiparticle residue Z and the 
mass renormalization are only weakly renormalized and 
remain finite at 5 C . These results suggest that the Lan- 



dau Fermi liquid parameters are strongly renormalized. 
In particular the instability is associated with Fq — > — 1 
as S is reduced toward 5 C , signaling the onset of phase- 
separation and/or superconductivity. 

We begin with the slave boson representation of the 
Hubbard model. In the infinite-U limit, the model de- 
scribes electrons with nearest neighbor hopping, t, on a 
2D square lattice, subject to the constraint that double 
occupancy on each site is prohibited. It is convenient to 
describe the projected Hilbert space in terms of a neutral 
spin-carrying fermion, fj a , creating the singly occupied 
site and a spinless charge-e boson, bi, keeping track of the 
empty site ||. The electron creation operator becomes 
c\ a = fj a bi. In the SU(N) generalization, the occupancy 
constraint thus translates into fi a fi<r + b}bi = N/2, where 
sum over repeated a — 1, . . . N index is implied. The 
partition function in the coherent state path integral for- 
mulation is 



Z = J VtfVbVpVfVXe fo L{T)dT , 
where the Lagrangian is given by 
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+ Y J ^Ulha + b^ l -N/2). 
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Here \i is a static Lagrange-multiplier enforcing the local 
constraint and fi is the chemical potential fixing an av- 
erage of S holes or n particles per site, i.e. < j\ a ji a >— 
N(l -S)/2 = n. The Lagrangian in Eq. (g) has a U(l) 
gauge symmetry, it is invariant under local U(l) transfor- 
mations: h — > bie %6 \ f ia — ► f lcr e' 9i , and A, — » A, - d T 8i. 
We choose the radial gauge § where the boson fields 
(bi,b\) are replaced by a real amplitude field r, while 
Xi is promoted to a dynamical field Aj(r). In this gauge, 
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the fermionic excitations can be identified with the Fermi 
liquid quasiparticles. 

To enable an 1 /N-expansion to the next-to-leading or- 
der, we write the boson fields in terms of static mean-field 
and dynamic fluctuating parts, 

n{T)=b[l + 5n{T)}, iXi(T)=X + i5Xi(T). (3) 

In the first part of the paper, we shall calculate b, A, 
together with the chemical potential p to the next-to- 
lcading order. Using these results, we then analyze 
the single-particle Green's function, determine the wave 
function renormalization Z and the quasiparticle mass 
renormalization and the compressibility. 

Substituting Eq. (jj|) into Eq. (j^), integrating out the 
fermions and the boson fields (<5r, <5A) to quadratic order 
in Eq. (|l|), we obtain the free energy F — —kThiZ to 
next-to-leading order in l/N, 
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where u) n is a fermion Matsubara frequency, ty. = 
— —fj-lk + A — p, with 7^. = cosk x + cosk y , and i*bos 
is the contribution due to boson fluctuations. The latter 
can be written in terms of the determinant of the inverse 
boson propagator matrix D , 



Fhaa = ^ In Pet J l (q,iv n ), 



(5) 



where v n is a boson Matsubara frequency. Note that 
in order to properly regularize the theory in the ra- 
dial gauge, DetZ? -1 should be evaluated on a discretized 
imaginary time mesh before taking the continuum limit 
in r The opposite sequence of operations will lead 

to unphysical ultraviolet singularities. We find, 

DetD^ 1 (q,iu n ) = P\\{q,iv n )P r r{<l,iv n ) - P 2 r (q,iv n ) 

+ 2b 2 [\ - e b (0)}S- P Xr (q, iv n )/iv n . (6) 



Here S~ = 



is a regularization fac- 



tor and eb(q) = A — 2tJ2lk-qnf(ek) with n/(e) the 
Fermi distribution function. P a p — N(H a p + B a p) are 
the fermion polarizations given by B rr — 2b 2 ei,(q)/N, 
B Xr = B rX = 2b 2 /N, Baa = 0, and 
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a/3 



E 



n/(£fc+) ~ n/(efc -W,^(M), (7) 



k £fc + 



where k± = fc ± q/2 and A = [—2tb 2 /N(^k + + 7fc_), i] are 
the boson-fermion vertices. 

The values of the parameters (b, A, p) are determined 
by minimizing the free energy in Eq. (^), leading to three 
self-consistent equations: 



OF _ () dF _ () dF 
db d\ ' dp 



(8) 



Solving these equations to leading order in l/N, where 
only the fermion contribution enters Eq. (Q), one recov- 
ers the results of Kotliar and Liu ||, namely, a bo- 
son condensate b 2 = b\ = N5/2 and a chemical po- 
tential shift A = Ao = 2t 7fc n /( £ fc)- This corre- 
sponds to a Fermi liquid phase with a quasiparticle dis- 
persion el = —{2wqI 'N)~fk + Ao — po an d a quasiparticle 
residue Z = b\ = N5/2 = m/m*. The compressibility 
k = dn/dp = Np/(l + 4tp\eo\), where p = J2k$( e k) an d 
peo = — Ylk 7fc^( e fc)- I* diverges as 5 — > 0, together with 
Z — > and m* — > oo, giving rise to a Brinkman-Rice jlj] 
metal-insulator transition. 
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FIG. 1. Phase shift A(q,v) at q = (2tt/3, 2tt/3) for 
S = 0.05, 0.15 and comparison to holon contributions. 

The effects of interactions between the quasiparti- 
cles and the incoherent Mott-Hubbard excitations enter 
through Fb os in Eq. (|4|) at the next-to-leading order in 
l/N HH . It is instructive to rewrite F bos in Eq. (j|) 
by converting the boson Matsubara sum into a contour 
integral distorted along the real axis, 



Fh 



-E 



dvA(q, v)n b (v), 



(9) 



where rif, is the Bose distribution function and 
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arctan 



ImDetL» _1 (q, v) 



KeT)etD- 1 {q,v) 



(10) 



can be considered as a many-body phase shift due to scat- 
tering of the fermions by particle-hole excitations. We 
have numerically calculated the phase shift A at T = 
from Eqs. ([l0|) and (^). Its general behavior is shown in 
Fig. 1 for a fixed wave vector q = (2n/3, 2ir/3) as a func- 
tion of frequency at different dopings. From intermediate 
to high frequencies, the scattering is in the unitary limit 
with A = 7T, indicating the existence of a collective mode 
which is pulled out of the particle-hole continuum at low 
frequency where A drops from ir to zero. Indeed, we find 
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that Det-D -1 has a branch-cut along the real axis cor- 
responding to the particle-hole continuum, and isolated 
poles corresponding to a collective mode which is well 
described by, 

lu 2 ~ c 2 [sm 2 (q x /2) + sin 2 (g,/2)] + e 2 (q), (11) 

where c oc St is the zero sound velocity and Cb{q) coincides 
with the original slave-boson dispersion. This mode has 
been identified as the "holon" in the t-J model ^|||- At 
small doping, the holon contribution, with uj* ~ ±£f,(g), 
dominates the particle- hole scattering as seen in Fig. 1. It 
disperses over the entire lower Hubbard band and carries 



where E k = -{2tb 2 /N) lk , S rX = e^ ' , and S\\ = 
(e~ lVn ° + e w ™° )/2 are regularization factors for D r \ 
and D\\ respectively. Without them, the theory in the 
radial gauge would be singular in the ultraviolet because 
D r \ and D\\ approach constants at large frequencies Q . 
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FIG. 2. Quasiparticle residue Z kp as a function of 
doping 8 in the TM direction. Inset: Slave boson conden- 
sate amplitude b/bo as a function of doping. 

Next we present the results of our numerical evalua- 
tions of Eqs.( |l2||L4| ), which were done on a 2D mesh of 
up to 60 x 60 points in the first quadrant of the Brillouin- 



the incoherent Mott-Hubbard spectral weight. Remark- 
ably, the holon contribution leads to a density-density 
correlation function in excellent agreement with that ob- 
tained from exact diagonalization of the t-J model on 
small clusters JlO[ ] - 

Now we solve the self-consistent equations in (|J) to 
next-to-leading order in 1/N including Fbos- Writing 
b = 60 + 61, A = Ao + Ai, and /i = /xq + fii, we find, 

^ = X 1 + ^H n (k,e k )S(e k ) + ^^, (13) 
P k N 



(14) 



(15) 

I 

zone using the micro-zone method and a frequency grid 
size as small as Aoj/t = 6/ 20 to ensure convergence. The 
result for the slave-boson condensate to next-to-leading 
order in 1/N is shown in the inset of Fig. 2 for N = 2. 
Interestingly, b vanishes at a doping 5* ~ 0.12. If we 
approximate the D rr in Eq. ( |l2| ) by the single holon 
mode in Eq. ([ll]) at small doping, we find an analyt- 
ical estimate 6/60 = 1 — 1/4-/V5, which vanishes at a 
5* = 1/4N = 0.125, in good agreement with the nu- 
merical result. 

It is important to note that at this order, the bo- 
son condensate is not simply related to the quasiparticle 
residue, in contrast to the case in the large-N limit. To 
determine the Fermi liquid coherence factor Z , we follow 
Refs. and write done the 1/N-resummation of the 
single-electron Green's function, 

ntvi \ b 2 [l + Z a (k,iLo n )} 2 ^ . 
G[k,iu n ) = . — - +b T,i(k,iu n ), (16) 

IbJn - £fc - E„(fc, lOJn) 

where E„ is given in Eq. (p"5|), Y, a and Uj are the anoma- 
lous part due to the boson condensate, and the incoherent 
part of the self-energies respectively. The latter are given 
by, to leading order in 1/N, 

E,-(fe, ilo) = -T*^2 G*o(fc + q, ico + w n )D rr {q,iv n ) (17) 
Y, a (k,iuj) = -T 22 G {k + q,iu + iv n )[D\ r (q,iv n )S r \ 
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]\j 2tb 2 1 

Ai = - ^ X] G o(k,iuj n )[£ n (k,iuj n ) - -^--^^ k ^ q D rr (q,w n )e~ lUn ° ] 

2t 2 

(k, iu! n ) + 2t\e \ X S„(fc, e k )8(e k ) + — 2J [D r \(q, iv n ) - D r \(q, 00)] . 

k,uj„ k P q,v n 

Here Gq 1 = iu> n — e k and E„(fc, iu) n ) is the usual self-energy to leading order in 1/N 

T, n (k,iu) n ) = X! lk-<jD rr (q, ii> n )e~ lUn0 - — ^ G (k + q,iu n + iVn)[D\\(q,iis n )S\\ 

q,iu n k,u n 

+2D\ r (q,iv n )S r \(E k + E k+q ) + D rr (q,iv n )(E k + E k+q ) 2 ], 
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-{E k + E k+q )D rr (q, w n )] . 



(18) 



The quasiparticle residue on the interacting Fermi sur- 
face can be obtained from Eq. (16), 



& 2 [l + Re£ a (fc F ,0)] 2 
[1 - <9Re£„(fc F , u)/du)\ u = ] ' 



(19) 



Thus Zk F can be finite even if b 2 is vanishing, provided 
that the reduction of the condensate is compensated by 
the contributions from the self-energies. Remarkably, 
this turns out to be the route followed by the 1/N- 
expansion. Expanding Eq. (|n]) to next-to-leading order 
in 1/N, one has 



A/N 



b 2 + 26 2 £ a (k F , 0) + b 2 dH n {k F , w)/&4 

A/N 



(20) 



In Fig. 2, Z£" is plotted as a function of doping in the 
TM direction. The 1 /N-corrections are clearly small and 
Zk F stays close to the large-N limit value. Within the sin- 
gle holon mode (Eq. ^) approximation, we found that 
the 1 / ^-correction to b 2 in Eq. ( ^p| ) is canceled out by the 
contributions from the self-energy terms, leaving Z^J^ 
weakly renormalized near 5*. Thus we conclude that 
while the boson condensate vanishes at 6*, the Fermi 
liquid coherence remains finite. 




FIG. 3. Electron chemical potential and the compress- 
ibility (inset) as a function of doping. 

We next turn to the compressibility of the model. In 
Fig. 3, the electron chemical potential fi = po + p\ is 
shown as a function of doping, which is strongly mod- 
ified from the N = oo result. The corresponding com- 
pressibility K = —d5/dfj, is shown in the inset of Fig. 3. 
At moderate dopings, n is approximately constant, but 
becomes strongly doping dependent as 5 is reduced. In- 
terestingly, there exists a critical doping, S c = 0.07±0.01, 
at which k diverges and becomes negative for 5 < 5 C . 
Thus, the Fermi liquid state becomes unstable below <5 C , 
while no singularity is present in Z kp To further un- 
derstand the nature of the instability, we have stud- 
ied the quasiparticle mass renormalization defined by 



m*/m = N*(Q)/N(0), where N(0) = p and AT*(0) are 
the bare (N = oo) and the renormalized (next-to-leading 
order in 1/N) quasiparticle density of states respectively. 
The numerical calculations of N* (0) show that while m* 
is enhanced in the doping range 0.05 < S < 0.2, it does 
not exhibit any singular behavior. A well-behaved N*(0), 
together with the general Fermi- liquid result, 



dn 
dp 



N*(0) 

TTfP' 



(21) 



suggests that the divergence of n is a result of the Lan- 
dau Fermi liquid parameter F® — > —1 at S c , indicative 
of phase-separation and/or superconducting instability 
p| . Note that the phase separation in the infmite-U 
case has a different origin than in models with strong 
antifcrromagnetic correlations. For one-hole, the ground 
state is known rigorously to be a Nagaoka state |14j of a 
saturated ferromagnet. For a finite density of holes, one 
expects ferromagnetic correlations to compete with the 
kinetic energy and whether the Nagaoka state remains 
stable is a question of great interest. Both numerical p"5| ] 
and analytical Jl6| results have shown that the uniform 
Nagaoka ferromagnetic state is unstable for any finite 
hole concentration. Our results naturally suggest a novel 
possibility that at low-doping the system phase separates 
into hole-poor ferromagnetic and hole-rich paramagnetic 
regions. In the presence of long-range Coulomb repulsion, 
we expect the p-wave pairing instability || enhanced 
by the tendency towards phase-separation to dominate 
p2],^3[ . We conclude that the breakdown of the Fermi 
liquid in our case is not due to a gradual reduction of 
the Fermi liquid coherence, but rather the enhanced in- 
teractions between the quasiparticles. This is the kind of 
Fermi liquid instability originally envisioned by Landau. 
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